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Abstract. This paper presents new results allowing an unknown non-Gaussian posi- 
tive matrix-valued random field to be identified through a stochastic elliptic boundary 
value problem, solving a statistical inverse problem. A new general class of non-Gaussian 
positive-definite matrix-valued random fields, adapted to the statistical inverse problems in 
C^ I high stochastic dimension for their experimental identification, is introduced and its prop- 

erties are analyzed. A minimal parametrization of discretized random fields belonging to 
this general class is proposed. Using this parametrization of the general class, a complete 
identification procedure is proposed. New results of the mathematical and numerical ana- 
►^ ' lyzes of the parameterized stochastic elliptic boundary value problem are presented. The 

^SJ ' numerical solution of this parametric stochastic problem provides an explicit approxima- 

c7^ ' tion of the application that maps the parameterized general class of random fields to the 

Oi . corresponding set of random solutions. This approximation can be used during the identi- 

fication procedure in order to avoid the solution of multiple forward stochastic problems. 
\l . Since the proposed general class of random fields possibly contains random fields which 

o 

approximation methods are introduced. In order to obtain an algorithm for constructing 
the approximation of a very high-dimensional map, complexity reduction methods are 
introduced and are based on the use of low-rank approximation methods that exploit the 
^^ ' tensor structure of the solution which results from the parametrization of the general class 

$H . of random fields. 

Notations 

A lower case letter, j/, is a real deterministic variable. 

A boldface lower case letter, y = (j/i, . . . , yjv) is a real deterministic vector. 

An upper case letter, F, is a real random variable. 

A boldface upper case letter, Y = (Yi, . . . , Ym) is a real random vector. 

A lower (or an upper) case letter between brackets, [a] (or [A\), is a real deterministic 

matrix. 

A boldface upper case letter between brackets, [A], is a real random matrix. 

E: Mathematical expectation. 
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Mjv,m(]K): set of all the (A'^ x m) real matrices. 
Mm(R): set of all the square {m x m) real matrices. 
M^(R): set of all the symmetric (n x n) real matrices. 

): set of all the symmetric positive-definite (n x n) real matrices. 
): set of all the skew-symmetric (m x m) real matrices, 
set of all the upper triangular (n x n) real matrices with strictly positive diagonal. 
0(7V): set of all the orthogonal [N x N) real matrices. 
Vm(R^): compact Stiefel manifold of [N x m) orthogonal real matrices. 
tr: Trace of a matrix. 

< y , z >2'. Euclidean inner product of y with z in M". 
||y||2: Euclidean norm of a vector y in R". 

\\A\\2: Operator norm subordinated to R": ||A|J2 = supnyn^^j ||Ay||2. 
Il^lli?: Frobenius norm of a matrix such that || A|||. = tr{[A]^ [A]}. 
[Im,n\- matrix in Mm,„(R) such that [/m,«]y = ^ij- 
[/,„]: identity matrix in M,„(R). 



1 Introduction 

The experimental identification of an unknown non-Gaussian positive matrix-valued ran- 
dom field, using partial and limited experimental data for an observation vector related to 
the random solution of a stochastic elliptic boundary value problem, stays a challenging 
problem which has not received yet a complete solution, in particular for high stochastic 
dimension. An example of a challenging application concerns the identification of random 
fields that model at a mesoscale the elastic properties of complex heterogeneous materials 
that can not be described at the level of their microscopic constituents (e.g. biological 
materials such as the cortical bone) , given some indirect observations measured on a col- 
lection of materials samples (measurements of their elastic response e.g. through image 
analysis or other displacement sensors). Even if many results have already been obtained, 
additional developments concerning the stochastic representations of such random fields 
and additional mathematical and numerical analyzes of the associated stochastic elliptic 
boundary value problem must yet be produced. This is the objective of the present paper. 

Concerning the representation of random fields adapted to the statistical inverse prob- 
lems for their experimental identification, one interesting type of representation for any 
non-Gaussian second-order random field is based on the use of the polynomial chaos 
expansion [67l [8] for which an efficient construction has been proposed in [25l [26l [17] , 
consisting in coupling a Karhunen-Loeve expansion (allowing a statistical reduction to be 
done) with a polynomial chaos expansion of the reduced model. This type of construction 
has been extended for an arbitrary probability measure [5S1 [351 [ISl [Ml [SB [20] and for 
random coefficients of the expansion [S5] . 

Concerning the identification of random fields by solving stochastic inverse problems, 
works can be found such as [15l [34l [35l [61] and some methods and formulations have 
been proposed for the experimental identification of non-Gaussian random fields |16| 
[27l [T2] [69] [T3l l4T] l60l [2] in low stochastic dimension. More recently, a more advanced 
methodology [57l [58] has been proposed for the identification of non-Gaussian positive- 
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definite matrix- valued random fields in high stochastic dimension for the case for which 
only partial and limited experimental data are available. 

Concerning the mathematical analysis of stochastic elliptic boundary value problems 
and the associated numerical aspects for approximating the random solutions (the for- 
ward problem) , many works have been devoted to the simple case of uniformly elliptic 
and uniformly bounded operators [Illllllll[ll[Slllllinil31].For these problems, ex- 
istence and uniqueness results can be directly obtained using Lax-Milgram theorem and 
Galerkin approximation methods in classical approximation spaces can be used. More 
recently, some works have addressed the mathematical and numerical analyzes of some 
classes of non uniformly elliptic operators [24l [28l |45l [10] , including the case of log-normal 
random fields. 

For the numerical solution of stochastic boundary value problems, classical non adapted 
choices of stochastic approximation spaces lead to very high-dimensional representations 
of the random solutions, e.g. when using classical (possibly piecewise) polynomial spaces 
in the random variables. For addressing this complexity issue, several complexity re- 
duction methods have been recently proposed, such as model reduction methods for 
stochastic and parametric partial differential equations [47j |48l [7] , sparse approximation 
methods for high-dimensional approximations |63l |46l [11] [6j [51] , or low-rank tensor ap- 
proximation methods that exploit the tensor structure of stochastic approximation spaces 
and also allow the representation of high-dimensional functions [181 l36l [50l [9l [lU [22l [43] . 

In this paper, we present a set of new results concerning: 

(i) the construction of a general class of non-Gaussian positive-definite matrix-valued 
random fields which is adapted to the identification in high stochastic dimension, 
presented in Section [51 

(ii) the parametrization of the discretized random fields belonging to the general class 
and an adapted identification strategy for this class, presented in Section [3] 

(iii) the mathematical analysis of the parameterized stochastic elliptic boundary value 
problems whose random coefficients belong to the parameterized general class 
of random fields, and the introduction and analysis of dedicated approximation 
methods, presented in Section [4] 

The numerical solution of the parameterized stochastic boundary value problems provides 
an explicit approximation of the application that maps the parameterized general class 
of random fields to the corresponding set of random solutions. This explicit map can be 
efficiently used in the identification procedure in order to avoid the solution of multiple 
stochastic forward problems. The general class of random fields possibly contain random 
fields which are not uniformly bounded, which requires a particular mathematical analysis 
and the introduction of dedicated approximation methods. Also, since the solution of 
this problem requires the approximation of a very high-dimensional map, we introduce 
complexity reduction methods using low-rank approximation methods that exploit the 
tensor structure of the solution which results from the particular parametrization of the 
general class of random fields. Note that this kind of approach for stochastic inverse 
problems has been recently analyzed in j52j in a particular Bayesian setting with another 
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type of representation of random fields, and with sparse approximation methods for the 
approximation of the high-dimensional map. 



2 General class of non-Gaussian positive-definite matrix-valued random 

fields 

2.1 Elements for the construction of a general class of random fields 

Let d ^ 1 and n ^ 1 be two integers. Let Z? be a bounded open domain of M'', let dD be 
its boundary and let Z) = Z? U dD be its closure. 

Let C+ be the set of all the random fields [K] = {[K(x)],x e D}, defined on a 
probability space {Q,T,V), with values inM+(]R). A random field [K] in C+ corresponds 
to the coefficients of a stochastic elliptic operator and must be identified from data, by 
solving a statistical inverse boundary value problem. The objective of this section is to 
construct a general representation of [K] which allows its identification to be performed 
using experimental data and solving a statistical inverse problem based on the explicit 
construction of the solution of the forward (direct) stochastic boundary value problem. 
For that, we need to introduce some hypotheses for random field [K] which requires the 
introduction of a subset of C+ . 

In order to normalize random field [K], we introduce a function x i-> [it](x)] from D 
into M+(M) such that, for all x in Z? and for all z in R", 

fcoll^ll^ s^ <[^(x)]z,z>2^ n-'^%\\z\\l, (2.1) 

in which k^ and ki are positive real constants, independent of x, such that < fcg < 
fcj < +00. The right inequality means that ||^(x)||2 ^ n^^^'^k^ and using ||Z^(x)|Jf ^ 
V^l|i^(x)||2 yields \\K{x)\\f s^ hi- Since sC [Z^(x)]jj and [Ki^)]'jj ^ \\K{^)\\f, it can 
be deduced that 

tr[^(x)] < fci, (2.2) 

in which fc j^ is a positive real constant, independent of x, such that A: i = nki. 

We then introduce the following representation of the lower-bounded random field [K] 
in C^ such that, for all x in Z3, 

[K(x)] = Y^[i(x)]^ {e[In] + [Ko(x)]} [L(x)] , (2.3) 

in which e > is any fixed positive real number, where [Z„] is the (n x n) identity ma- 
trix and where [L{x.)] is the upper triangular (n x n) real matrix such that, for all x in 
D, [K{x)] = [L{^)f [i(x)] and where [Kq] = {[Ko(x)],x G D} is any random field in C+ . 

For instance, if function [K] is chosen as the mean function of random field [K], that is 
to say, if for all x in D, [Z^(x)] = Z;{[K(x)]}, then Eq. ([231) shows that Z;{[Ko(x)]} must 
be equal to [/«]) what shows that random field [Kq] is normalized. 

Lemma 1 If [Kq] is any random field in C^ , then the random field [K], defined by 
Eq. i2.3\) . is such that 
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(i) for all X in D, 

||K(x)||f ^-^(V^e+||Ko(x)|lf) a.s. (2.4) 

I + e 

(ii) for all z in M" and for all x in D, 

fcj|z||2 <^ <[K(x)]z,z >2 a.s, (2.5) 

in which k^ — kQe/{l + e) is a positive constant independent of x.. 
(iii) for all x in D, 

||[K(x)]-i||j.^ ^^^ + ^\ r[X(x)]-i a.s, (2.6) 

which shows that, for all integer p ^ 1, {[K(x)]^^,x € D} is a p-order random 
field, i.e., for all x in D, _B{||[K(x)]^^||'^} < +oo and in particular, is a second- 
order random field. 

Proof (i) Taking the Frobcnius norm of the two members of Eq. (|2.3p . using |lL(x)||ir 
= ||[L(x)]^||f = ^ytI[K_{-x.)] and taking into account Eq. ([2?2l) . we obtain Eq. (pT4)) . 



(ii) Eq. p.Sp can easily be proven using Eq. (|2.3p and the left inequality in Eq. ()2.ip . 
(iii) We have ||[K(x)]-i||f ^ |l[L(x)]-i||f ||[L(x)]-^||^||[K,(x)]-i||j. jn which [K,(x)] = 
(e[ J„] + [Ko(x)])/(l + e). Since [Ko(x)] is positive definite almost surely, for x fixed in 
D, we can write [Ko(x)] = [$(x)] [A(x)] [$(x)]-^ in which [A(x)] is the diagonal random 
matrix of the positive- valued random eigenvalues Ai(x), . . . , A„(x) of [Ko(x)] and [*(x)] 
is the orthogonal real random matrix made up of the associated random eigenvectors. 
It can then be deduced that [Ke(x)]-i = (1 + e) [*(x)] (e[/„] + [A(x)])-i [*(x)]^ and 
consequently, for x in D, ||[K,(x)]-i||2, = tr{[K,(x)]-2} ^ (l + e)^ ^;Li(e + A,(x))-2 < 
(l + £)2n/£2. Since |1[L(x)]-1||f ||[i(x)]-^||F = tr[X(x)]-i, we deduce Eq. I^Bl- □ 

Lemma 2 // [Kq] is any random field in C^ , such that, for all x in D, 

(i) ||Ko(x)||i? ^ /3o < oo almost surely, in which /3o is a positive-valued random vari- 
able independent of:x., then the random field [K], defined by Eq. 112.3]) . is such that, 
for all Ti. in D, 

|1K(x)|1f ^ /3<+oo a.s, (2.7) 

in which /3 is the positive-valued random variable independent of x such that (3 = 
k^{V^e + l3o)/{l + e). 
(ii) i?{||Ko(x)|||,} < +00, then [K] is a second-order random field, 

E{\\K{x)\\l} < +00 . (2.8) 

Proof Eqs. (PT)) and ([^ are directly deduced from Eq. (pTi)) . D 



Remark 2.1 If /3o is a second-order random variable, then Eq. (j2.7p implies Eq. ()2.8p . 
However, if Eq. (|2.8p holds for all x in D, then this equation does not imply the existence 
of a positive random variable /3 such that Eq. (|2.7p is verified and a fortiori, even if /3 
existed, this random variable would not be, in general, a second-order random variable. 
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2.2 Representations of random field [Kq] as transformations of a 
non- Gaussian second-order symmetric matrix- valued random field [G] 

Let C^''^ be the set of all the second-order random fields [G] = {[G(x)],x g £)}, defined 
on probability space (8,T, P), with values in Mfj(M). Consequently, for all x in D, we 
have£;{||G(x)|||,}<+oo. 

In this section, we propose two representations of random field [Kq] belonging to C^ , 
yielding the definition of two different subsets of C^: 

• Exponential type representation. The first type representation is written as [Ko(x)] = 
exp]^([G(x)]) with [G] in C^' and where expm denotes the exponential of symmetric 
square real matrices. It should be noted that random field [G] is not assumed to 
be Gaussian. If [G] were a Gaussian random field, then [Kg] would be a log-normal 
matrix-valued random field. 

• Square type representation. The second type representation is written as [Ko(x)] = 
[L(x)]"^ [L(x)] in which [L(x)] is an upper triangular (n x n) real random matrix, for 
which the diagonal terms are positive-valued random variables, and which is written 
as [L(x)] = [£([G(x)])] with [G] in C^-^ and where [G] ^ [C{[G])] is a well defined 
deterministic mapping from M^(R) into the set of all the upper triangular [n x n) real 
deterministic matrices. Again, random field [G] is not assumed to be Gaussian. If for 
all 1 ^ i ^ i' ^ n^ [d-Jjj' are independent copies of a Gaussian random field and for 
a particular definition of mapping [£([0])], then [Kq] would be the set SFG^ of the 
Non-Gaussian matrix- valued random field previously introduced in j55| . 

These two representations are general enough, but the mathematical properties of each 
one will be slightly different and the computational aspects will be different. 

2.2.1 Exponential type representation oj random field [Kq] 

It should be noted that for all symmetric real matrix [A\ in M^(M), [B] — expi^([A]) is 
a well defined matrix belonging to M^(R). All matrix [A\ in M^(M) can be written as 
[A\ — [$] [/i] [^Y' ill which [/i] is the diagonal matrix of the real eigenvalues fii, . . . ,fin 
of [A] and [$] is the orthogonal real matrix made up of the associated eigenvectors. We 
then have [B] = [$] expp^([/x])[$]-^ in which expj^ ([//]) is the diagonal matrix in M+(R) 
such that [expM[([/x])]jfe = e'"' Sjk ■ 

For all random field [G] belonging to C^'^ , the random field [Kg] — expfy[([G]) defined, 
for all X in _D, by 

[Ko(x)] = expM([G(x)]) , (2.9) 

belongs to C^ . If [Kg] is any random field given in C^ , then there exists a unique random 
field [G] with values in the M^(M) such that, for all x in Z?, 

[G(x)]=logM([Ko(x)]). (2.10) 

in which log]y[ is the reciprocity mapping of expj^ which is defined on M+ (M) with values 
in M|(M), but in general, this random field [G] is not a second-order random field and 
therefore, is not in C^'^ . The following lemma shows that, if a random field [Kq] in C^ , 
satisfies additional properties, then there exists [G] in C^'^ such that [Kq] = expiyi([G]). 
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Lemma 3 Let [Ko] be a random, field belonging to C^ such that for all x Cz D, 

ii;{||Ko(x)||2,}<+^ , i?{||[Ko(x)]-i|||,}<+^. (2.11) 

Then there exists [G] belonging to C"^^ such that [Kq] = expj^([G]). 

Proof For x fixed in D, we use the spectral representation of [Ko(x)] introduced in 
the proof of Lemma [T] The hypotheses introduced in Lemma [3] can be rewritten as 
ii;{||Ko(x)|||} = E;=i E{K,(xf} < +00 and E{\\ [Ko(x)]-i|||} = E;=i E{A,{x)-^} < 
+C)o. Similarly, it can easily be proven that i?{||G(x)|||n} = X]'i=i E{{logAj{x})'^}. Since 
n is finite, the proof of the lemma will be complete if we prove that, for all j, we have 
E{{logAj{x))'^} < +CX3 knowing that E{Aj{x.y^} < +oo and E{Aj{x)'^} < +oo. For j 
and X fixed, let P{dX) be the probability distribution of the random variable Aj(x). We 
have E{{\ogAj{x)f} = /^ (log A)^ P{dX) + /+°°(log A)^ P{dX). For < A < 1, we have 
(logA)^ < A~^, and for 1 < A < +oo, we have (logA)^ < A^. It can then be deduced 
that £;{(log Aj(x))2} < J^ A-2 P(dA) + /+°° A^ P{dX) < /+°° A-^ p(dA) + /+°° A^ P{dX) 
= E{A,{x)-^} + ^{Aj(x)2} < +00. D 

Remark 2.2 The converse of Lemma |3] does not hold. If [G] is any random field in C^'^ , 
in general, the random field [Kq] = exp[^([G]) is not a second-order random field. 

Proposition 1 Let [G] be a random field belonging to C^' and let [Kg] be the ran- 
dom field belonging to C+ such that, for all x in D, [Ko(x)] = exp]^([G(x)]). Then the 
following results hold: 

(i) For all x in D, 

||Ko(x)|lf ^ V^ell^f'')"- a.s, (2.12) 

i?{||Ko(x)|||,}^nii;{e2||G(x)||.}^ (2.13) 

i?{||logM[Ko(x)]|i|,}<+(X3. (2.14) 

(ii) // ||G(x)||i? ^ /3g < +00 almost surely, in which Pq is a positive random variable 
independent of ii., then ||Ko(x)||i? =^ /3o < +00 almost surely, in which /3o is the 
positive random variable, independent of x, such that j3q — ^/n e^'^ . 

Proof For (i), we have || exp[^(G(x))||2 ^ ell*^''''!!^. Since [Ko(x)] = expiy(G(x)) and 
since ||Ko(x)||f ^ V^I|Ko(x)||2, it can then be deduced that ||Ko(x)|li^ ^ v^el'*^'")!!^ 
and then i?{||Ko(x)|||,} ^ ni?{e^ll*^('''ll^}. Since [G] is a second-order random field, the 
inequality (|2.14p is deduced from Eq. (|2.10p . For (ii), the proof is directly deduced from 
the inequality ((^1^ . D 



2.2.2 Square type representation of random field [Kq] 

Let g 1—^ h{g;a) be a given function from M in IR+, depending on one positive real 
parameter a. For all fixed a, it is assumed that: 
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(i) h{.]a) is a strictly monotonically increasing function on M, which means that 
h{g] a) < h{g'] a) if — cxd < g < g' < +00; 

(ii) there are real numbers < c/j < +00 and < c^ < +cxd, such that, for all g in M, 
we have h[g; a) ^ Ca+ Chg^ ■ 

It should be noted that Ca ^ /i(0, a). In addition, from (i) it can be deduced that, for 
g < 0, we have h{g;a) < h{0,a) ^ Ca < Ca + cug'^. Therefore, the inequality given in 
(ii), which is true for g < 0, allows the behavior of the increasing function g i— > h{g]a) 
to be controlled for g > 0. Finally, the introduced hypotheses imply that, for all a > 0, 
g 1-^ h{g; a) is a one-to-one mapping from R onto R^ and consequently, the reciprocity 
mapping, v 1— )■ h~^{v; a), is a strictly monotonically increasing function from M+ onto R. 

The square type representation of random field [Kg] belonging to C+ is then defined as 
follows. For all x in D, 

[Ko(x)]=L([G(x)]), (2.15) 

in which {[G(x)],x g £)} is a random field belonging to C^"^ and where [G] 1-^ h{[G]) is 
a measurable mapping from M^(R) into M+(R) which is defined as follows. The matrix 

[Ko] = L([G]) e M+(R) is written as 

[Ko] = [Lf [L] , (2.16) 

in which [L] is an upper triangular (n x n) real matrix with positive diagonal, which is 
written as 

[L]=C{[G]), (2.17) 

where [G] t-^ C{[G]) is the measurable mapping from M^(M) into Mj/(R) defined by 

mG])]n' ^ [G],r ,l^j<f^n, (2.18) 



mG])]n = \/H[G]n,a,) , 1 ^ j ^ n , (2.19) 

in which ai, . . . , a„ are positive real numbers. 

If [Kg] is any random field given in C^, then there exists a unique random field [G] 
with values in the Mfj(R) such that, for all x in 13, 

[G(x)]=L-i([Ko(x)]), (2.20) 

in which L"-"^ is the reciprocity function of L, from M+(R) into Mfj(R), which is explicitly 
defined as follows. For all 1 ^ j ^ j' ^ n, 

[G(x)],y =[/:-i([L(x)])],y , [G(x)],., = [G(x)],y . (2.21) 

in which [L] 1— )> £~^([L]) is the unique reciprocity mapping of C (due to the existence of 
V I— >■ /i~^(i'; a)) defined on Mj;'(R), and where [L(x)] follows from the Cholesky factoriza- 
tion of random matrix [Ko(x)] = [L(x)]^ [L(x)] (see Eqs. ^^U^i and ((TTB)l ). 



Example of function h 

Let us give an example which shows that there exists at least one such a construc- 
tion. For instance, we can choose h = /i^'''^, in which the function h^^'^ is defined 
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in [55] as follows. Let be s = 6/\/n + 1 in which (5 is a parameter such that < 
S < i/(n + l)/(n — 1) and which allows the statistical fluctuations level to be con- 
trolled. Let be qj = 1/(2 s^) + (1 - j)/2 > and h'^^ig-^a) = 2 s^ Ff]^ (Fw (g / s)) with 
Fwiw) — J^^ — i= exp(— iw^) dw and F-^^{u) = 7 the reciprocal function such that 

FrAl) = u with Fr„(7) = /J' rR ^""^ e"* di and T{a) ^ /+°° f-^ e"* di. Then, for 
all j = 1, . . . ,n, it can be proven [55] that g 1— >■ h'^^'^{g;aj) is a strictly monotonically 
increasing function from M into K"*" and there are positive real numbers c^ and Ca ■ such 
that, for all g in R, we have h'^^'^{g; aj) ^ Ca + c/i g^. In addition, it can easily be seen 
that the reciprocity function is written as h'^^'^~ {v;a) — sF^^^(Fr^(u/(2s^)). 

Proposition 2 Let [G] be a random field belonging to Cj{^ and let [Kq] be the random 
field belonging to C^ such that, for all x in D, [Ko{x)] = L([G(x)]). We then have the 
following results: 

(i) There exist two real numbers < 70 < +00 and < 71 < +cx) such that, for all x 
in D, we have 

||Ko(x)|lf ^7o+7i|lG(x)|l| a.s, (2.22) 

E{\\Ko{^)\\f}<+^, (2.23) 

//^{||G(x)|||,}<+oo, then E{\\Ko{x.)\\l} < +00 . (2.24) 



(ii) //||G(x)||i? ^ /3g < +00 a.s, in which Pq is a positive random variable independent 
of:ic, we then have |jKo(x)||i? ^ /3o < +00 almost surely, in which /3o is the positive 



random variable, independent of x, such that /3o = 7o + 7i 0c 



Proof (i) Since [Ko(x)] = [L(x)]^ [L(x)], wehave |1Ko(x)||f =$ ||L(x)|||, with ||L(x)|||, = 
E,/^([G(x)],,;a,) + E,</[GWLy- Therefore, ||Ko(x)||f «; E, Ca, + max{c^, 1/2} 
^ , [G(x)]^j, which yields inequality (I2.22p with 70 = X^i "^qj ^^^^ 71 — max{c;i, 1/2}. 
Taking the mathematical expectation yields i*^{ II Ko(x) II i?} ^ 7o + 7ii?{||G(x)|||.} < +00 
because [G] is a second-order random field. Taking the square and then the mathematical 
expectation of Eq. (|2:22)) yields Eq. ((2:24)) . 
(ii) The proof is directly deduced from Eq. (|2.22p . D 



2.3 Representation of any random field [G] in C^ "^ 

As previously explained, we are interested in the identification of the random field [K] 
which belongs to C+, by solving a statistical inverse problem related to a stochas- 
tic boundary value problem. Such an identification is carried out using the proposed 
representation of [K] (defined by Eq. (|2.3p ) as a function of the random field [Kg] 
which is written either as [Ko(x)] — expj^([G(x)]) (see Eq. ()2.9p ) or as [Ko(x)] = 
[£([G(x)])]^ [/:([G(x)])] (see Eqs. (|2l6l) and (f2l7l) ). In these two representations, [G] is 
any random field in C^''^, which has to be identified (instead of [K] in C+). Consequently, 
we have to construct a representation of any random field [G] in C^^^ . Since any [G] in 
C^''^ is a second-order random field, a general representation of [G], adapted to its iden- 
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tification, is the polynomial chaos expansion for which the coefficients of the expansion 
constitutes a family of functions from D into M^(R). 

It is well known that a direct identification of such a family of matrix- valued functions 
cannot easily be done. An adapted representation must be introduced consisting in choos- 
ing a deterministic vector basis, then representing [G] on this deterministic vector basis 
and finally, performing a polynomial chaos expansion of the random coefficients on this 
deterministic vector basis. The representation is generally in high stochastic dimension. 
Consequently, the identification of a very large number of coefficients must be done. It 
is then interesting to use a statistical reduction and thus to choose, for the deterministic 
vector basis, the Karhunen-Loeve vector basis. Such a vector basis is constituted of the 
family of the eigenfunctions of the compact covariance operator of random field [G] . 



2.3.1 Covariance operator of random field [G] and eigenvalue problem 

Let [Go(x)] = ii^{[G(x)]} be the mean function of [G], defined on D with values in 
M^(E). The covariance function of random field [G] is the function (x, x') i— )■ Cg(x, x'), 
defined on Z? x Z?, with values in the space M^(R) ® M|(]R) of fourth-order tensors, such 
that 

Cg(x,x') = i?{([G(x)] - [Go(x)]) ® ([G(x')] - [Go(x')])} . (2.25) 

It is assumed that [Go] is a uniformly bounded function on £>, i.e. 

||Go||oo := CSS supxg£,||Go(x)||F < +oo (2.26) 

and that function Cq, is square integrable on D x D. Consequently, the covariance oper- 
ator Covg, defined by the kernel Gq, is a Hilbert-Schmidt, symmetric, positive operator 
in the Hilbert space L^(D,M^(M)) equipped with the inner product, 

« [G,] , [G,] »= / tr{[G,(x)]^ [G,(x)]}cix, (2.27) 

Jd 

and the associated norm |||[Gi]||| =^ [Gi] , [Gt] ':$>^^^. The eigenvalue problem related to 
the covariance operator Covg, consists in finding the family {[Gi(x)],x G D}i^i of the 
normalized eigenfunctions with values in M^(M) and the associated positive eigenvalues 
{ciji^i with o"! J^ (72 ^ . . . — > and X]i=T '^i "^ +oo, such that 

/ Gg(x,x') : [G,(x')] dx' - a, [G,(x)] , i^l, (2.28) 

Jd 

in which {Gg(x,x') : [Gi(x')]}fe£ = J2k'e'{CG{^,x')}kek'e' {[Gi(x')] }*,/£/. The normalized 
family {[Gi(x)],x e D}i^i is a Hilbertian basis of L^(I?,M^(R)) and consequently, 
<C [Gi] , [Gj] 3>= Sij. It can easily be verified that 

+00 



S{||G(x)||2,}^||Go(x)||2,+^a,||G,(x)||2,<+oo, Vx e Z? , (2.29) 

i=l 

/ Z;{|lG(x)|l|,}dx-|||Go||P+Va,<+oo. (2.30) 

Jd ,_i 



i=l 
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We now give some properties which will be useful later. Since Cq is a covariance func- 
tion, we have ||Cg(x,x')|P s^ tr{CG(x,x)}x tr{CG(x',x')} in which ||Cg(x,x')||2 = 
Y.kik'i'{'^G(.^^^')}kek'i' and tr{CG(x,x)} = Y.ke{^Gi^^^)}keke ^ 0. We can then de- 
duce the following Lemma. 

Lemma 4 // x i— > tr{CG(x, x)} is integrable on D, then Cg is square integrable on 
D X D. If:x.i-^ tr{CG(x, x)} is bounded on D, then Cg is bounded on D x D (and thus 
square integrable on D x D) and the eigenfunctions x t-^ [Gi(x)] are bounded functions 
from D into M„(R). For i ^ 1, we then have HGiHoo — £ss sup^^j^\\Gi{x.)\\p < +oo. 

Using additional assumptions about the kernel Cg , more results concerning the decrease 
rate of eigenvalues {(Tiji^i can be obtained. 

Lemma 5 We have the following results for d = 1 and for d ^ 2. 

(a) For d = 1, it is proven (see Theorem 9.1 of Chapter II of fTOJ /) that, if for a given 
integer /i ^ 1 and for all integer a such that 1 ^ a ^ /x, the functions (a;, a;') i— > 
d^Ccix^ x')/dx''^ are bounded functions on D x D, then X^iJT '^i "^ ^"°° 
for any (^ > 0. 

(b) For d ^ 2 (finite integer), a useful result is given by Theorem 4 of fSTl. Let us 
assume that D has a sufficiently smooth boundary dD. For x in D, let h be such 
that X + h and x + 2h belong to D . We then define the difference operator Ah 
such that, for all x' fixed in D, AhCG(x, x') = Cg(x -1- h, x') — Cg(x, x') and A^ 
is defined as the second iterate of Ah. For all x' fixed in D, for a given positive 
integer fi such that /i ^ 1 and for a given real C such that < C ^ 1 > let be 

||Cg(-,x')||, = 5: sup||P«Cg(x,x')||+ E sup^i^^^i^, 

in which X'°Cg(x,x') ^ al°ICG(x,x')/ax"i ■■■dx'^" with \a\ ^ ai + . . . + aa- 
Therefore, if Cg is continuous on D x D such that sup^/g^ ||Cg(-, x')||^ < +oo, 
we then have the following decrease of eigenvalues: ai — 0(1^^^^^^^'^'''^). 

(c) For d ^ 2 (finite integer), a similar result to (b) in Lemma [5| can be found in 
f44) / under hypotheses weaker than those introduced above, but in practice, it seems 
much more difficult to verify these hypotheses for a given kernel Cg ■ 



Comment about the eigenvalue problem 

Any symmetric (n x n) real matrix [C] can be represented by a real vector w of dimension 
riw = n{n -{- l)/2 such that [C] = ^(w) in which Q is the one-to-one linear mapping 
from R"™ in Mfj(M), such that, for 1 ^ i ^ j ^ n, one has [G]ij — [C]ji — Wk with 
k = i + j{j — l)/2. Let £^ be the set of all the R"™-valued second-order random fields 
{W(x),x e D} defined on probability space {<d,T,'P). Consequently, any random field 
[G] in C"^'^ can be written as 

[G]=g(W), W = g-\[G]), (2.31) 
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in which W is a random field in C\ . The eigenvalue problem defined by Eq. (|2.28p can 
be rewritten as 

/ [Cw(x,x')] wXx')dx' = f7,w*(x) , i ^ 1, (2.32) 

JD 

in which [Cw(x,x')] = £;{(W(x)-wO(x))(W(x')-wO(x'))'^} and where wO = g-i([Go]). 
For 1^1, the eigenfunctions [Gi] are then given by [Gi\ = 5(w'). 



2.3.2 Chaos representation of random field [G\ 

Under the hypotheses introduced in Section [2.3.1[ random field [G] admits the following 
Karhunen-Loeve decomposition, 

[G(x)] - [Go(x)] + ^ ^^- [G,(x)] r;, , (2.33) 

1=1 

in which {r]i}i^i are uncorrelated random variables with zero mean and unit variance. 
The second-order random variables {rji]i^i are represented using the following polyno- 
mial chaos expansion 

-t-oo 

i=i 
in which {S^jfegN is a countable set of independent normalized Gaussian random variables 
and where {^j}j>i is the polynomial chaos basis composed of normalized multivariate 
Hermite polynomials such that E{'^ j{{'E.k}k£n)'^ i'{{'^k}k£n)} = 6jj'- Since E{rjir]ii} = 
dii> , it can be deduced that 

+ C30 

Y^yiyi>=6u'. (2.35) 



2.4 Random upper bound for random field [K] 

Lemma 6 If^^J^ ^/^\\Gi\\oo < -l-oo, then for all x in D, 

\\G{yi)\\ F ^ fie < +00 a.s., (2.36) 

in which /3q is the second-order positive-valued random variable, 

+ 00 

|3G = \\Go\\oo+Y.V^^\\G^\\ooh\: ^{/3|} < +00 . (2.37) 

i=l 

Proof The expression of /3g defined in Eq. (|2.37l) is directly deduced from Eq. (I2.33p . 
The random variable Pa can be written as /3g = ||Go||oo + P- Clearly, if E{P^} = 
E,.»'v^v^l|G'Joo||G,'||ooS{|?7*l \Vi'\} < +00, then £;{/32;} < +00. Wehave£;{|7?,| \t],,\} ^ 
^/EWiVEWi = 1 and thus E0^} ^ (ES V^I|G.||oo)^. D 
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Remark 2.3 It should be noted that, if HGiHoo < Coo < +00 for all i ^ 1 with Coo 

■\+oo 



independent of i, we then have X^iJT v^ll^'iljoo < +00 if: 

(i) for d = 1, the hypothesis of (a) in Lemma [5] holds for /i = 2. The proof is the 
following. We have J2t=^ ^/a7||Gi||oo < Coo J2t=i V^- Since CTj -> for i ^- +00, 
there exists an integer io ^ 1 such that, for all i ^ io, we have ai < 1. For ^ — 2, 
then there exits < C < 1/10 such that J^t^o V^ ^ ^t^o ^i ^ ^°°' ^^^'^^ 

yields X^iJT V^ < +00, and consequently, X]i=^ ^oo y/o^ < +00. 

(ii) for d ^ 2 (finite integer), the hypothesis of (b) in Lemma [5] holds for jjl = d. The 
proof is then the following. We have X^iJT v^ll^'ll°o ^ '^oaYlit^ V^ ^'^'^ ^'^^ 
^i = d, ^ = 0(r i"'^/(2'^)) and consequently, for < C sJ 1, ES V^ < +°°- 

The previous results allow the following proposition to be proven. 

Proposition 3 Let cr^ and [Gi] be defined in Section [2.3.1\ If^^^ y^WGiWoo < +00, 
then for all x in D, 

||K(x)||i. s; /3< +00 a.s, (2.38) 

in which j3 is a positive-valued random variable independent of x. Let j3q be the second- 
order positive-valued random variable defined by Eq. {2.31^ . 

(i) (Exponential type representation) // [K(x)] is represented by Eq. 112. 3\) with Eq 



[K(x)] = Y^[i(x)]^ {£[/„] + expM([G(x)])} [L(x)] , (2.39) 

then, 13 = fciVn(e + e'^'^)/(l + e). 

(ii) (Square type representation) // [K(x)] is represented by Eq. 112. 3\) with Eqs. 112.16]) 
and I^J^ , 

[K(x)] = Y^[^(x)]^ {£[/„] + [C{[Gi^)])f [/:([G(x)])]} [L(x)] , (2.40) 

then, (3 = ki{-\/n£ + 70 + 71 0q)/{1 + e) in which 70 and 71 are two positive and 
finite real numbers. In addition the random variable 13 is such that 

E{I3}< +00. (2.41) 

Proof Proposition [3] results from Lemma [2J Propositions [T] and [5J and Lemma IHl D 

2.5 Approximation of random field [G] 

Taking into account Eqs. (|2.33p to (12.351) . we introduce the approximation [G*^'"-^^] of 
the random field [G] such that, 

in 

[G(™^^) (x)] = [Go (x)] + ^ V^ [G. (x)] Ty, , (2.42) 

i=l 

N 

^. = Ey'*^(^)' (2.43) 
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in which the {'i'jIjLi only depends on a random vector H = (Si, . . . , ^Ng) of A''^ inde- 
pendent normalized Gaussian random variables Si, ... , Sat^ defined on probability space 
(B,T, P). The coefficients yl are supposed to verify J2j=iyi Vi' — ^n' which ensures 
that the random variables, {?7i}™i, are uncorrelated centered random variables with 
unit variance, which means that E{riiriii} = Sw. The relation between the coefficients 
can be rewritten as 

[yf[y] = [lm], (2.44) 

in which [y] E Mjv,m(K) is such that [y]ji ~ yj ioi 1 ^ i ^ m and 1 ^ j ^ N. Introducing 
the random vectors rj — (771, . . . ,77™) and *(H) = (^i(H), . . . , *Ar(H)), Eq. (|2.43p can 
be rewritten as 

r, = [y]^*(H). (2.45) 

Equation p.44p means that [y] belongs to the compact Stiefel manifold 

V,„(R^) ^{[y]e MAr^,„(R) ; [yf [y] = [/,„] } . (2.46) 

With the above hypotheses, the covariance operators of random fields [G] and [G'-™'^^] 
coincide on the subspace spanned by the finite family {[Gi]}™ 1. 



3 Parametrization of discretized random fields in the general class and 

identification strategy 

Let us consider the approximation {[G''"'^-'(x)],x G D} of {[G(x)],x e D} defined by 
Eqs. (|2.42p to (|2.44p . The corresponding approximation {[K('"'^^(x)],x G D} of random 
field {[K(x)],x e D} defined by Eq. (|239)) or by Eq. (I2;i0)) . is rewritten, for aU x in £>, 
as 

[K(™.^)(x)]=/C(™^^)(x,H,[y]), (3.1) 

in which (x, y, [y]) M- ]Qbn,N) ^^^ y^ j^j-j jg ^^ mapping defined on _D x R^s x Vm(M^) with 
values in M+(M). 

• The first objective of this section is to constructed a parameterized general class of 
random fields, {[K'-™'^^(x)],x G -D}, in introducing a minimal parametrization of the 
compact Stiefel manifold Vm(R^) with an algorithm of complexity 0{Nm^). 

• The second objective will be the presentation of an identification strategy of an optimal 
random field {[K(™'^)(x)],x S D} in the parameterized general class, using partial 
and limited experimental data. 



3.1 Minimal parametrization of tiie compact Stiefel manifold V„i(R'^) 

Here, we introduce a particular minimal parametrization of the compact Stiefel manifold 
V,„(R^) using matrix exponentials (see e.g. [T]). The dimension of V„i(R^) being v = 
mN — m{m+l)/2, the parametrization consists in introducing a surjective mapping from 
R'' onto Vm(R^)- The construction is as follows. 

(i) Let [a] be given in Vm(R^) and let [a±] G MAr.jv_m(R) be the orthogonal com- 
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plement of [a], which is such that [a aj_] is in 0{N). Consequently, we have, 

[a]'^ [a] = [/„i], [aj_]'^ [aj_] = [/jv-m], [a±]'^ [a] = [OAr_,„,,„] . (3.2) 

The columns of matrix [a±] can be chosen as the vectors of the orthonormal basis 
of the nuh space of [a]'^ . In practice [5U], [a±] can be constructed using the QR 
factorization of matrix [a] = [Qn] [R-m] and [a±] is then made up of the columns 
j = m + l,...,N oi[QN] eO{N). 

(ii) Let [A] be a skew-symmetric (m x m) real matrix which then depends on m{m — 
l)/2 parameters denoted by zi, . . . , Zrn(m-i)/2 and such that, ioi I ^ i < j ^ m, 
one has [A\ij — —[A\ji — Zk with k = i + (j — l)(j — 2)/2 and for 1 ^ i ^ m, 
[A]u = 0. 

(iii) Let [B] be a {{N''m)xm) real matrix which then depends on {N—m)m parameters 



denoted by 2;„(„_i)/2+l: ■ • ■ : ^m(m-l)/2+(Af-m)m 



and such that, for 1 ^ z ^ iV — m 

with k = i + {j - 1){N - m). 



and 1 ^ j' ^ m, one has [B\i 

Introducing z — (zi, . . . , z^) in R"^, there is a one-to-one linear mapping S from 

W into Mf,f (R) X MN-m.rr 



R) such that 

{[A],[i?]} = 5(z), (3.3) 

in which matrices [A\ and [B] are defined in (i) and (ii) above as a function of z. 

Let t > be a parameter which is assumed to be fixed. Then, for an arbitrary [a] in 
Vm(R^), a first minimal parametrization of the compact Stiefel manifold Vm(R^) can 
be defined by the mapping Mwa from R"^ onto Vm(R^) such that 



[y] = MiaA'Z') ■■= [a a_L] {expM(i 



A -B^ 
B 



)}[l 



N,m\ 



(3.4) 



In the context of the present development, we are interested in the case for which N ^ m 
and possibly, in the case for which N ^ m with N very large. The evaluation of the 
mapping M[a] defined by Eq. p.4p has a complexity 0{N^). 

We then propose to use a second form of minimal parametrization of the compact 
Stiefel with a reduced computational complexity. This parametrization, which is derived 
from the results presented in [19], is defined by the mapping M[a] from R"^ onto Vm(R^) 
such that 



[y] = -^[a](z) := [a Q] {expM(t 



A 
R 



^R^ 




)}[h 



(3.5) 



in which [a Q] G Mat 2m (R). The matrix [Q] is in Vm(R ) and [R] is an upper triangular 
(m X m) real matrix. These two matrices are constructed using the QR factorization of 
the matrix [aA_] [B] e M7v,m(R), 

[a^] [B] = [Q] [R] . (3.6) 

The evaluation of the mapping Mu] defined by Eq. (|3.5p has a complexity 0{Nm^). 
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3.2 Parameterized general class of random fields and parameterized random 

upper bound 

Using Eq. p.Sp . the parameterized general class of random field {[K'^™'^''(x)], x e D} is 
then defined as 

/C :- {/C(™'^)(-,H,X[„](z)) ; z e M'^} . (3.7) 

It should be noted that, for z = 0, [y] = A4[a]{0) = [a], which corresponds to the random 
field [K*^'"'^^] = /C(™^^)(-, H, [a]). The following proposition corresponds to Proposition|3] 
for the approximation [K^™'^''] of random field [K]. 

Proposition 4 The random field [K^™'^)] = IC'-"^-'^\-,S,,M[a]{z)), z e R" , is such that 

IIK^-'^'lIf s^7(H,z)<+oo (3.8) 

almost surely and for all z e R"^, where 7 : R^« x W — > K is a measurable positive 
function. For the square type representation of random fields, there exists a constant 7, 
independent on N and z, such that 

£;{7(H,z)} s;7< +00 forallzeMy. (3.9) 

Moreover, */ X^iJT v^ll^*lloo ^ +0O; wii/i (7; and [Gi] defined in Section \2.3.1\ p.9p is 
satisfied for a constant^ independent of m. 

Proof Following the proof of Lemma IH we obtain 

■m 

||G("^'^)(x)||f < IIGolloo +^ V^||G.||ooh^(H,z)| := 5(H,z), (3.10) 

with 77 = (771, . . . ,7]m) = A^[o](z)'^*(H). Using Proposition [31 we then obtain Eq. (|3.8p 
with 7 = fciA/n(e + e'')/(l+e) for the exponential type representation and 7 = ki{^/ne + 
7o + 71 '5^)/(l + £) for the square type representation. Using E{Tfi} = 1, it can be shown 
that EiS"^} ^ 2||Go||L + 2 (E"i %/^l|G»||oo)^ := Cm, where C,„ < +00 is independent 
on N and z. Therefore, for the square type representation, we obtain Eq. p.9p with 
7 = fci(V"e + 70 + 71 Cm)/(1 + £)• Moreover, if Y.'t^ ^i\\Gi\\ 00 < +00, then C,„ s^ 
Coo < +00 andEq. ([31]) holds for 7 = ^^(Vrie + 7o + 71 Coo)/(l + e)- □ 

3.3 Brief description of the identification procedure 

Let B be the nonlinear mapping which, for any given random field [K] = {[K(x)], x G D} 
introduced and studied in Section [21 associates a unique random observation vector 
■jjobs _ g([K]) with values in R'^ots, The nonlinear mapping B is constructed in solving 
the elliptic stochastic boundary value problem as explained in Section [4] We are then 
interested in identifying the random field [K] using partial and limited experimental data 
set u*^"?'^, . . . , u™P''^'=''P in R^obs (^niobs and I'exp are small). In high stochastic dimension 
(which is the assumption of the present paper), such a statistical inverse problem is an 
ill-posed problem if no additional available information is introduced. As explained in 
|57 [ 158 ) . this difficulty can be circumvented (1) in introducing an algebraic prior model 
(APM), [K'^''^], of random field [K], which contains additional information and satisfying 
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the required mathematical properties and (2) in using an adapted identification proce- 
dure. Below, we summarize and we adapt this procedure to the new representations and 
their parametrizations of random field [K] that we propose in this paper. The steps of 
the identification procedure are the following: 

Step 1 

Introduction of a family { [K^''^(x; w)] , x e D} of prior algebraic models for random field 
[K]. This family depends on an unknown parameter w (for instance, w can be made up 
of the mean function, spatial correlation lengths, dispersion parameters controlling the 
statistical fluctuations, parameters controlling the shape of the tensor- valued correlation 
function, parameters controlling the symmetry class, etc). For fixed w, the probability 
law and the generator of independent realizations of the APM arc known. For examples 
of algebraic prior models of non-Gaussian positive-definite matrix-valued random fields, 
we refer the reader to [SSj for the anisotropic class, to l62j for the isotropic class, to [30] 
for bounded random fields in the anisotropic class, to [21] for random fields with any 
symmetry class (isotropic, cubic, transversal isotropic, tetragonal, trigonal, orthotropic), 
and finally, to |32| for a very general class of bounded random fields with any symmetry 
properties from the isotropic class to the anisotropic class. 

Step 2 

Identification of an optimal value w°p' of parameter w using the experimental data set, 
the family of stochastic solutions U°'"*(w) = B([K*™(-; w)]) and a statistical inverse 
method such as the moment method, the least-square method or the maximum likelihood 
method [53l [59l iMl [57]). The optimal algebraic prior model {[KO^™(x)],x e D} := 
{[K'^'''^(x; w°P')],x G D} is then obtained. Using the generator of realizations of the 
APM, i/KL independent realizations [X^^)], . . . , [ivT^'^L)] of random field [K°*™] ^^^ ^^ 
generated with Vgj^ as large as it is desired without inducing a significant computational 
cost. 

step 3 

Choice of a type of representation for random field [K] and, using Eq. (12. 3p with Eq. (|2.10p 
or with Eq. (I2.20[) . the optimal algebraic prior model {[G°*''^(x)], x e D} of random field 
[G] is deduced. For all x e D, [G°^™(x)] = logM[([K°A™(x)]) for the exponential type 
representation and [G°^™(x)] = L^^([Kq*''^(x)]) for the square type representation, 
with 

[Kr^(x)] = (1 + £)[i(x)]-^ [K°^™(x)] [L(x)] - e[/„] . (3.11) 

It is assumed that random field [G°'^''^] belongs to C^'^ ■ From the ly^L independent real- 
izations [-f^^^^], . . . , [/'C'^'^kl)] of random field [K°^''^], it can be deduce the Vyll independent 
realizations [G^^)], . . . , [G^'^kl)] of random field [G° 



-"OAPMl 



Step 4 

Use of the i/kl independent realizations [G^^-'], . . . , [G^'^kl)j of random field [G°'^''^] and 
use of the adapted statistical estimators for estimating the mean function [Gq*'"'^] and 
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the covariance function Cgoapm of random field [G°'^''^] (see Section 12.3. ip . Then, calcu- 
lation of the first m eigenvalues cti ^ ... ^ Om and the corresponding eigenfunctions 
[Gi], . . . , [Gm] of the covariance operator Covgoapm defined by the kernel Cgoapm. For a 
given convergence tolerance with respect to m, use of Eq. (|2.42p to construct independent 
realizations of the random vector ryO^P^ = [rfl'^^^, . . . , ry^f '^) such that 

m 

[QOAPM(m)(^)j = [Gg*™(x)] +J2V^, [G,(x)] ?7°^™ . (3.12) 

i=l 

For i — 1, . . . ,m, the z/kl independent realizations t]^ , . . . , 77^*''' of the random variable 
j^oAPM ^j.g calculated by 

r/f ^ = ^ « [G(^)] - [GSn JG,] » , ^-1,...,;.KL. (3.13) 



Step 5 

For a given convergence tolerance with respect to N and A^g in the polynomial chaos 
expansion defined by Eq. ()2.45p . use of the methodology based on the maximum likelihood 
and the corresponding algorithms presented in [57] for estimating a value [yo] E V„i(K^) 
of [y] such that 7?°'^''^ = [yof *(H). 

Let us examine the following particular case for which the algebraic prior model [K'^''^] 
of random field [K] is defined by Eq. ^^ with either Eq. ([2l0|) or Eq. ^^, in which 
j-QAPMj jg Qjjosen as a second-order Gaussian random field indexed by D with values 
in M„(R). Therefore, the components rji, . . . ,rifn of the random vector t] defined by 
Eq. (|2.45p are independent real-valued normalized Gaussian random variables. We then 
have Ng = m. Let us assume that, for 1 ^ j ^ m ^ N , the indices j of the polynomial 
chaos are ordered such that ^j(H) = Ej. It can then be deduced that [yo] E Vm(M^) is 
such that [yolji = 5ij for 1 ^ i ^ ?ti and 1 ^ j ^ N. 

Step 6 

With the maximum likelihood method, estimation of a value z of z G W^ for the parame- 
terized general class /C defined by Eqs. (|3.ip and (|3.7p . using the family of stochastic solu- 
tions U°''^(z) = S(/C("^^) {■,a,M[y^] (z))) and the experimental data set u'^"?-!, . . . , u'=''P'''-p. 
Then, calculation of [y] — M[yg]{z). 

Step 7 

Construction of a posterior model for random field [K] using the Bayesian method. In 
such a framework, the coefficients [y] of the polynomial chaos expansion rj = [y]'^ ^(■3) 
(see Eq. (|2.45p ') are modeled by a random matrix [Y] (see [SB]) as proposed in [SH] and 
consequently, z is modeled by a R'^-valued random variable Z. For the prior model [YP"""^] 
of [Y], here we propose 

jYPrior] ^ j^ ^_^ (2P™'') , ZP"" - centered Gaussian vector , (3.14) 

which guaranties that [YP"""^] is a random matrix with values in Vm(R^) whose statis- 
tical fluctuations are centered around [y] (the maximum likelihood estimator of the set 
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of coefficients [y] computed at step 6). The Bayesian update allows the posterior dis- 
tribution of random vector ZP"'*' to be estimated using the stochastic solution U"*'^ = 
B(/C('"'^n-, S, X[^](ZP"™))) and the experimental data set u'=''P'\ . . . , u'=''P'''-p. 

Finally, it should be noted that once the probability distribution of ZP"'"*' has been 
estimated by Step 7, Vkl independent realizations can be calculated for the random field 
[GP°^'(x)] = [G°^™(x)] +E"iV^[G«(x)]r/P'''' in which rjP°'' = [YP"**']^ *(H) and 
where [YP°^'] = A^[g](ZP°''*'). The identification procedure can then be restarted from 
Step 4 replacing [G^'^''^] by [GP"'^]. 



4 Solution of the stochastic elliptic boundary value problem 

Let D C M'' be a bounded open domain with smooth boundary. The following elliptic 
stochastic partial differential equation is considered, 

- div([K] • VC/) = / a.e. inD, (4.1) 

with homogeneous Dirichlet boundary conditions (for the sake of simplicity) . The random 
field {[K(x)],x e D} belongs to the parameterized general class of random fields. 



/C = {/C("'^)(.,H,A^[,](z));zeR''}, 



introduced in Section 13.21 This stochastic boundary value problem has to be solved 
at steps 6 and 7 of the identification procedure described in Section 13.31 respectively 
considering z as a deterministic or a random parameter. We introduce the map 

u: Dx R^s X M'" ^ R 

such that U — m(x, H,z) is the solution of the stochastic boundary value problem for 
[K(x)] = /C^'"'^^(x, H,7V{[a](z)). The aim is here to construct an explicit approximation 
of the map u for its efficient use in the identification procedure. 

In this section, a suitable functional framework will first be introduced for the definition 
of the map u. Then, numerical methods based on Galerkin projections will be analyzed 
for the approximation of this map. Different numerical approaches will be introduced de- 
pending on the type of representation of random fields (square type or exponential type) 
and depending on the properties of approximation spaces. For the two types of represen- 
tation of random fields (exponential type or square type). Finally, we will briefiy describe 
complexity reduction methods based on low-rank approximations that exploit the tensor 
structure of the high-dimensional map u and allows its approximate representation to 
be obtained in high dimension. That makes affordable the application of the identifica- 
tion procedure for high-dimensional germs H (high Ng) and high-order representation of 
random fields (high v). 

4.1 Analysis of the stochastic boundary value problem 

We denote by F = T^g (E) Fj, a product measure on R'^ := R^s x R"^, where FjVg is the 
probability measure of random variable H and where F^^ is a finite measure on R". Up to 
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a normalization, Fj, is considered as the probability measure of a random vector Z. We 
denote by [C] : D x M^" x M" ^ M+(R) the map defined by 

[C](x,y,z)=/C("^^)(x,y,A^[,](z)), 

and such that [K(x)] = [C](x, H, Z) is a ^(H, Z)-measurable random field. Sometimes, the 
random field {[C](x, H, Z), x e D} will be denoted by {[C(x)], x g D}. For a measurable 
function h : R^' x M'' — > M, the mathematical expectation of h is defined by 

Erih) = E{h{S, Z)} = / h{y, z) F(dy, dz) . 



Lemma 7 Under the hypotheses of Proposition[^ there exists a constant a and a positive 

measurable function 7 : R^s x M'' ^- R such that, for T-almost all (y,z) in R^s x R"^, 

we have 

n/ ^ -r -r < [C](x, y, z)h, h >2 , . „, 

< a ^ ess mi ml ——^ , 4.2) 

xsD heR"\{o} ||h||2 

< [C](x,y,z)h,h>2 ... ,, „x 

ess sup sup ., .,2 ^7(y,z)<oo. (4.jj 

xGD heK"\{0} l|h||2 

Moreover, for the class of random fields corresponding to the square type representation, 
we have Eri^j) ^ 7 < +00, that means 7 S Lp(R^). 

Proof Lemma [TJn) gives the existence of the lower bound a = k^. Proposition |4] yields 
||[C](x, y,z)||2 ^ ||[C](x,y, z)!!^ < 7(y,z) < 00, with 7 a measurable function defined 
on M^9 X R'^. For the square type representation of random fields, property p.9p implies 

Er{-f) = 4. i?r„^{7(H,z)}F,(dz) ^7. D 

We introduce the bihnear form C(-, sy^z) : Hq{D) x Hq{D) -^ R defined by 

Ciu,v;y,z) ^ f Vv [C](-,y,z)Vudx. (4.4) 

JD 

1/2 

Let introduce || • ||^i = (J^ |V(-)p dx) the norm on H^{D), and || • \\h-^ the norm on 
the continuous dual space H~^{D). 



Strong-stochastic solution 

Proposition 5 Assume f G H^^{D). Then, for T-almost all (y,z) in R^s x R"^, there 
exists a unique u{-,y,z) G Hq{D) such that 

C{u{-,y,z),v;y,z) ^ f{v) for all v e H^{D) , (4.5) 

and 

M;y,z)\\H.^-\\f\\H-^. (4.6) 

" a 
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Proof Lemma [7] ensures the continuity and coercivity of bilinear form C(-,-;y, z) for 
F-almost all (y, z) in R^s x M.'^. The proof follows from a direct application of the Lax- 
Milgram theorem. D 



Weak-stochastic solution 

Let beL^(R^) = L^^ {R^o)(SLUW) und X = H^{D) (g) Ll{R^'). Then X is a Hilbert 
space for the inner product norm || • ||x defined by 



J D 



We also introduce the spaces X'^^^'^'"^ {s,r £ N) of functions v : D x Ri^ ^ R with 
bounded norm 



Lemma 8 We have a^^'^\\v\\x ^ a|kllx(c) ^ Ikllxcc^) =^ ll^'llx(^c;) «$ ll^'llx(-'^)' "-^^ 
therefore 

x(y') c x('^f^) c x(c') c x(c) C X 

with dense embeddings. 

Proof The inequalities satisfied by the norms are easily deduced from the properties of 
[C] (Lemma [7]). Then, it can easily be proven that X'-'' ' is dense in X , which proves the 
density of other embeddings. D 

Let introduce the bilinear form a : X x X ^i' M. defined by 

a{u,v) = Er{C{u,v)) = f { j Vv [C](-,y,z)Vudx') r(dy,dz) , 

JRA" \Jd J 

and the linear form F belonging to the continuous dual space X' of X, defined by 

{F,v)^Er{f{v))= [ f{v{-,y,z))r{dy,dz). 
From Lemma [71 it can easily be deduced the 

Lemma 9 a : X x X ^>- W is a symmetric bilinear form such that: 
(i) a is continuous from X^'-^ ^ x X to M., 

\a{u, v)\ sC ||u||^(c2, ||i;|U V(w, v) e X(c') X X . (4.7) 

(ii) a is continuous from X'*^-* x X^'-^' to M, 

Hu,v)\^\\u\\xiaA\v\\xic^) y{u,v) e X^^^ X X^^K (4.8) 

(iii) a is coercive, 

a(w,t;) >a||f||3f Vw G X . (4.9) 
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Now we introduce a weak form of the parameterized stochastic boundary value problem: 

Find u & X such that 

(4.10) 

a[u,v)^F{v) yveX. 

We have the following result. 

Proposition 6 There exists a unique solution u Cz X to problem (j4.10p . and \\u\\x ^ 
i||F||x'- Moreover, u verifies (|4.5p T -almost surely. 

Proof The existence and uniqueness of solution can be deduced from a general result 
obtained in '45'. Here we provide a short proof for completeness sake. We denote by 
A : D{A) C X —¥ X' the linear operator defined by Au = a{u, •), D{A) being the 
domain of A. We introduce the adjoint operator A* : D{A*) C X —> X' defined by 
{u,A*v) = {Au,v) for ah u € D{A) and v € D{A*), with D{A*) = {v e X;3c > 
such that | {Au,v) \ < c||m|Jx for all u € D{A)}. The continuity property (|4.7p implies 
that A is continuous from X'''^ ' to X' and therefore, D{A) D X^*^ ^ is dense in X 
(using Lemma [8|). That means that A is densely defined. The coercivity property ()4.9p 
implies that ||Au|jx' ^ a||'y||x for all t; G X, which implies that A is injective and the 
range R{A) of A is closed. It also implies that A* is injective and R{A*) is closed, and 
by the Banach closed range theorem, we then have that A is surjective, which proves 
the existence of a unique solution. Then, using the coercivity of a, we simply obtain 
\\ur^^^a{u,u) = j^F{u)^i;\\F\\x>\\u\\x. D 

Proposition 7 The solution u of problem (|4.10|) is such that: 

(i) u e x^^\ 

(ii) uex^-^^^ if-feL}.{W'), 

(iii) uex^f^^ if-/e Ll{R^'). 



Proof Using Eq. (|4.5I) with v — u and Eq. (j4.6p . we obtain J^ Vu • [CjVudx = f{u) ^ 
ll/llff-ill'-'llffi =^ -ill/llH-i- Taking the expectation yields ||m||^(c) s^ ill/ll?/-i < oo, 
which proves (i). If 7 e Lf, ||w||^(,c) = Er{j J^Vu ■ [C]VMdx) s^ ^||/||^-i£;r(7) < 00, 
which proves (ii). Finally, if 7 € L^, we have ||u|l^(^2) = Erij'^WufHl) ^ i7ll/llH-i-Sr(7^) < 
00, which proves (iii). D 



4.2 Galerkin approximation 

Galerkin methods are introduced for the approximation of the solution u of problem 
(|4.10p . Let Xn C X be an approximation space such that Xn '■= Xq^p — Vq ® Wp 
with Vq C Hq{D) (e.g. a finite element approximation space) and Wp C Lp(R'^) (e.g. a 
polynomial chaos approximation space). The Galerkin approximation un G X^ of u is 
defined by 

a{uN,VN) = F{vn) ^vn e Xn . (4-11) 
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Note that the coercivity property of bihnear form a on Xn x Xj^ (Lemma [9jiii)) ensures 
the existence and uniqueness of a solution un to problem (|4.1ip . The convergence of 
Galerkin approximations is now analyzed in different situations corresponding to the 
different types of representation of random fields (exponential or square type), to different 
underlying measures F, and to different choices of approximation spaces. We first analyze 
the case where Xn C X^'^'> which leads us to a natural strategy for the definition of 
a convergent sequence of approximations. Then, we analyze the use of more general 
approximation spaces that do not necessarily verify Xn C X^'^\ 

4.2.1 Case Xn C X^'^) 

Proposition 8 Assuming X^ C X'^"'^ C X'^^\ the solution un G Xn of (14. lip verifies 

\W-un\\x(c) ^ inf \\u-vn\\x(c), (4.12) 

and \\u—un\\x^^= inf ||m — wjv||^(c) . (4.13) 

Therefore, if {XN}NeN C X^^> is a sequence of approximation spaces such that UngnXn 
is dense in X, then there exists a subsequence of Galerkin approximations {uN}NeN which 
converges to u in the X^^' -norm and in the X-norm. 

Proof Using the Galerkin orthogonality property of un and the continuity of a (Lemma 
ISTii)) yield 

II u- un\\x{c) = a{u- un,u- Un) = a{u -un,u~vn) ^ \\u - mjv||x(C) ||w - VN\\x(ch 

for all vn G Xn- Inequality (|4.12|) is obtained by taking the infimum over all vn G Xn- 
Then, inequality (|4.13p is obtained by using the coercivity of a (Lemma [DJiii)). D 

Let us now analyze the condition Xn C X^''^ of Proposition [H with Xn = Vg (X" Wp- 
Due to the tensor product structure of Xn = V^ (g) Wp, the condition Xn C X^'^^ is 
equivalent to ^/jWp C Lp(M^), that means Er{ip'^^) < +oo for all ip G Wp- 

Use of weighted approximation spaces (for both types of representation of 
random fields) 

A possible approximation strategy consists in choosing weighted approximation spaces 
Wp = {7~^/^(/3; iy9 G Wp} with Wp C L^{W)- For example, if the measure F has finite 
moments of any order, Wp can be chosen as a classical (piecewise) polynomial space with 
degree p, e.g. 

Wp = Pp(R^s) (E) FpiM.") = spanjy^z'^; a G N^«, /3 G N", |a| < p, \j3\ s: p}- 

Indeed, for aU if = ^j-'^l'^y'^^P g Wp, 



JRf 



+00, 



This approximation strategy is adapted to both types of representation of random fields 
(exponential or square type) since it does not require any assumption on 7. However, the 
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use of weighted polynomial approximation spaces Wp leads to non classical computational 
treatments. 



Use of classical approximation spaces (for square type representation of 
random fields) 

More classical approximation spaces can be used provided some conditions on 7. 

Lemma 10 // 7 G ip(R^) and Wp C Lp(]R^) with some r > 1 and s > 2 such that 
I + i = 1, then Xn C X('^). In particular Xn C X('^) i/7 G i^ (R'') and Wp C Lf{Rt'). 

Proof For r > 1, letting s' = | such that p- + ^ = 1, we have 

Eriip'j) ^ {Eriip'^')}'/^' {Er{Y)V^' = {Eri^nV^^ {ErhlV^^ < 00. 
The proof for r = 1 and s = 00 is straightforward. D 

Lemma[TO]allows us to analyze the approximation when using the square type represen- 
tation of random fields, for which 7 G Lp(M''). Indeed, in this case, LemmafTOl implies that 
Xn C X^''^ if Wp C L^(R^). In particular, this condition is satisfied if T has a bounded 
support and Wp is a (possibly piecewise) polynomial space, e.g. Wp = Pp(R^9)(g)Pp(M''). 
We note that F = F^v <8> Fi, has a bounded support if (i) the random germ H for the 
representation of random fields in the class K. = {/C''"'^-'(-, H, A^[(j](z)) ; z e M'^} has a 
bounded support supp{TN^), and (ii) the support of the measure Fi, on the parameter 
space R'' is bounded. This latter condition implies that the map u only provides the solu- 
tion to the boundary value problems associated to a subset of random fields in /C (those 
corresponding to parameters z £ supp(T^)). In other words, this first result shows that 
when using the square type representation of random fields with a germ H with bounded 
support, the approximation is possible using classical polynomial approximation spaces. 

Note that for having X^ C X^''\ weaker conditions on approximation spaces could 
be obtained by looking further at the properties of 7 and the measure F. In particular, 
if Wp is the tensor product of polynomial spaces with (total or partial) degree p, that 
means Wp = Pp(R^s) ^ Pp{W), then Xn C X^'^' if for all y^z'^ g Wp, 

Er{y^"z^Pj{y,z)) = / y^'-z^fj{y,z)r{dy,dz) < +00, (4.14) 

which is a condition on 7 and on measure F. The following result justifies the applicability 
of classical (piecewise) polynomial spaces when using the square type representation of 
random fields. 

Proposition 9 For the square type representation of random fields, if Wp is the tensor 
product of (possibly piecewise) polynomial spaces and ifT = Tn <8) F^ is such that prob- 
ability measures Fjv and T^^ admit moments of any order, then Xn = Vq (E) Wp C X^'^K 
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Proof Following the proof of Proposition HI we obtain 7 ^ ki{y/ne + 70+71(5^), with S 
defined by Eq. p.lOl) and such that 

m r, m 

i=l i=l 

with go = 2||Go||^ and g, = 2^™ 1 ^^IIGiH^. From Eq. ^M, we have 77 = M^*(y) 
where '^'(y) is a vector of polynomials in y and [j/]"^[y] — [Im]- The latter condition 
implies \yj\ ^ I for all 1 ^ i < to and 1 ^ j ^ TV. Therefore, rjf < iJ2]=i yl^j{y)Y ^ 

^E7=i *j(y)': and <5(y,z)2 <; go +5imiVEf=i *.(y)' := Q(y)> where Q(y) is a poly- 
nomial in y. If the measures F^r and F,^ have finite moments of any order, then any 
(piecewise) polynomial function on R^ is in ip(R^). Therefore, if Wp is a space of (piece- 
wise) polynomial functions on R'' , then for all ip G Wp, we have 

Er{ip{y, z)27(y, z)) < Er{ip{y, z)2(l + Q(y)) < -t-00, 

from which we deduce that Vq ® Wp C X'-^^ D 

Since measures with bounded support have finite moments of any order, Proposition 
[HI is consistent with the first conclusions of Lemma [TUl Moreover, we have that for the 
square type representation of random fields, if Z is chosen as a Gaussian random variable 
(that means F^ is a Gaussian measure) or a random variable with bounded support, then 
the classical Hermite polynomial chaos space associated with a Gaussian germ H can be 
used. Note that the use of a measure Fi, whose support is R"^ allows to explore the whole 
class of random fields /C with the single map u. 

Remark 4.1 The case of log- normal random fields corresponds to a particular case 
of the exponential type representation for which the random field [G] is Gaussian and 
represented using a degree one Hermite polynomial chaos expansion with a Gaussian germ 
H. In this case, it can be proven that if Wp is the tensor product of (possibly piecewise) 
polynomial spaces and if F,y admit moments of any order, then Xjv = Vq (S) Wp C X^'^' . 
This justifies the use of polynomial approximation spaces when considering the particular 
case of log-normal random fields. However, this result does not extend to other random 
fields with exponential type representation. 



4.2.2 Case of general approximation spaces Xn 

In order to handle both types of representation of random fields (exponential type and 
square type) with a measure F with a possibly unbounded support, we here propose 
and analyze an approximation strategy which consists in using truncated approximation 
spaces. We consider a family of approximation spaces {Vq}q^i such that UqVq is dense 
in Hq{D), and a family {Wp}p^i such that UpWp is dense in Lp(R^). Classical (possibly 
piecewise) polynomial spaces with degree p can be chosen for Wp. Let X^r be defined as 
Vq (8) Wp. Then, for r > 0, we introduce the approximation space 

Wp = {ij^rv; <p e VFp} , 
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where /^^^(y, z) = 1 if 7(y, z) ^ r and if 7(y, z) > r. It can be proven that Ur>o Up^i 
Wp is dense in Lp(]R^). Let X^ be defined as Vq (g) W^ . For ip — I-^^tV i^^ W^p j we 
have i;r(7V'^) ^ T£;r((/3^) < oo, that means Xl; := Vq W^ C X't) for all t > 0. The 
Galerkin approximation of u in XJf is denoted by u]^ . 

Proposition 10 The Galerkin approximation u'^ G XJ^ of u satisfies 

ll"-w]vllx(c) ^^ inf ||M-t;||^ + ||u/^>r||^(c), (4.15) 

vEXn 

and there exists a sequence of approximation spaces X'^,-, — Vg(T-) ® W^,s such that 

||u-u]v(^)||j(-(C) ^^ as T^oo. (4.16) 

Proof We have u G X^^'^ (Proposition [T]) and XJ^ C X^'') C X^*^). We first note that 
for all V € XJ^, 

\\{u-v)\W^c) = \\{u - v)Li^rfx(a) +||(w-u)/^>r|lx(C) 
S; ||(M-w)/7^r|lx(T) + l|w/7>r|lx(C) 
S; t\\{u - V)l^<^r\\]c + \Wll>r\\\(C)- 

Then, using Proposition [8] (Eq. ()4.13p ). we obtain 

ll"-"]vllx(C) =^ inf ll^-'^llxic) =^ 'T inf \\{u - v)I^^r\\\ + Wul^yrfxic) 
= T inf \\iu-v)I^^r\\x + hIi>T\\x(c) 
^r inf |Jm-u||| + ||m/^>^||^(c) 

vEXn 

Since ||w||x<c) is bounded, the second term converges to as r — )• oo. Also, provided that 
UisfXpf is dense in X, we can define a sequence of spaces ^Ar(7-) such that 

T inf ||m — w||^— >0 as r — > oo, 
which ends the proof. D 

4.3 High-dimensional approximation using low-rank approximations 

The approximation of the map u Cz X requires the introduction of adapted complex- 
ity reduction techniques. Indeed, when introducing approximation spaces Vq C H^(D) 
and Wp = Wy (g) W^ C L^^ (M^s) (g) L^_^(M''), the resulting approximation space 
Xn — Vq (E) Wp may have a very high dimension dim{Vq) x dim{Wp). For classical 
non adapted constructions of approximation spaces Wp, the dimension of Wp has typi- 
cally an exponential (or factorial) increase with Ng and v (e.g. with polynomial spaces 
W^ — Pp(M^9) and Wp — Pp(R'^)). Therefore, the use of standard approximation tech- 
niques would restrict the applicability of the identification procedure to a class of random 
fields with a germ H of small dimension (small Ng) and a low order in the representation 
of random fields (small v) . 
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Low-rank approximations 

Low-rank tensor approximation metliods can be used in order to reduce tlie complexity 
of the approximation of functions in tensor product spaces |33] . They consist in approxi- 
mating the solution in a subset A4 C X^ of low-rank tensors. Different tensor structures 
of space Xjv can be exploited, such as Xn = Vg (g) Wp or X^r = Vq (S) W^ (g) Wp. Also, if 
the measures T^^ and T^, are tensor product of measures, then W^ and Wp can be cho- 
sen as tensor product approximation spaces W^ = 'Sij^l^Wp'^ ' and W^ — ^k^iWp ' . 
Low-rank approximation methods then consist in approximating the solution in a subset 
A4 C X]\f of low-rank tensors, which is a subset of low dimension in the sense that M 
can be parameterized by a small number of parameters (small compared to the dimen- 
sion of Xn). Different low-rank formats can be used, such as canonical format. Tucker 
format, Hierarchical Tucker format or more general tree-based tensor formats (see e.g. 
[33]). These tensor subsets can be formally written as 

M ^{v^FM{wi,...,wey,wi eR''\...,weeW'} 

where F_m is a multilinear map with values in Xn and where the Wk S K'^'" {k = 1, . . . ,i) 
are the parameters. The dimension of such a parametrization is X]/c=i '''=• ^^ ^^ example, 
M can be chosen as the set of rank-m canonical tensors in Xn = Vq(S> W^ (S> Wp, defined 

hy M^{v^ j2Zi < '^ w'f ® <; < € Vq, wj e wy, < e w^}. 



Algorithms for the low-rank approximation of un € Xn 

The Galerkin approximation ujv G Xn of u is the unique minimizer of the strongly convex 
functional J : Xjv -> R defined by J{v) — ^a{v,v) — F{v). A low-rank approximation 
Ur (z A4r of UN Can then be obtained by solving the optimization problem 

ram J(v) — min J(Fm{wi, . . . ,we)). 

v£M wi£R''i,...,we£R''e 

This problem can be solved using alternated minimization algorithms or other optimiza- 
tion algorithms (see [33]). Also, greedy procedures using low-rank tensor subsets can 
be introduced in order to construct a sequence of approximations {u'^jfe with u''+^ = 
yfc _|_ ^fc+i ^^^ ^fe+i g j^ defined by J{u'' + w'=+^) = min^g^n J{u'' + v). We refer 
to [22] for the analysis of a larger class of greedy algorithms in the context of convex 
optimization problems in tensor spaces, and to |50) for the practical implementation of 
some algorithms in the context of stochastic parametric partial differential equations. 



4.4 Remarks about the identification procedure 



The identification procedure has been presented in Section 13.31 In view of the numerical 
analysis presented above, the following remarks can be done. Different approximations 
of the map u should be constructed at the different steps of the identification procedure. 
Indeed, the quality of approximations of u clearly depends on the choice of measure Fi, 
on the parameters space 'R'^. It is recalled that the parametrization of the class /C = 
^](^{m,N) j-.^ a,M.[a] (z)) ; z e M"^} of random fields depends on the set of parameters [a] = 
A^[a](0). The map u depends on the choice of [a] but it is clearly independent on the 
choice of measure T^, provided that the support supp{T^) = W^ . However, the Galerkin 
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approximation u^q being optimal with respect to a norm that depends on the measure 
Fj^, the quahty of the approximation clearly depends on measure Ty. For that reason, the 
choice of measure F^ can be updated throughout the identification procedure in order to 
improve the approximation of u as explained at the end of Section [3731 Step 6 and step 7 
described in Section 15751 can then be clarified as follows. 

At Step 6, the parametrization of /C around [a] = [yo] is used and an approximation of 
the corresponding map u in X is constructed using low-rank tensor methods. The use of 
this explicit map allows a fast maximization of the likelihood function to be performed, 
yielding an estimation z of z (or equivalently [y] of [y] ) . 

At Step 7, the map u constructed in Step 6 could be used for the subsequent Bayesian 
update. However, it is also possible to construct a new map that takes into account 
the result of the maximum likelihood estimation. Therefore, the parametrization of /C 
around [a] = [y] is used and the boundary value problem is solved again in order to 
obtain an approximation of the corresponding map m(-, •, •; [y]) in X . This map can then 
be efficiently used for solving the Bayesian update problem. 



5 Conclusion 

In this paper, we have presented new results allowing an unknown non-Gaussian positive 
matrix-valued random field to be identified through a stochastic elliptic boundary value 
problem, solving a statistical inverse problem. In order to propose a constructive and 
efficient methodology and analysis of this challenging problem in high stochastic dimen- 
sion, a new general class of non-Gaussian positive-definite matrix- valued random fields, 
adapted to the statistical inverse problems in high stochastic dimension for their exper- 
imental identification, has been introduced. For this class of random fields, two types 
of representation are proposed: the exponential type representation and the square type 
representation. Their properties have been analyzed. A parametrization of discretized 
random fields belonging to this general class has been proposed and analyzed for the two 
types of representation. Such parametrization has been constructed using a polynomial 
chaos expansion with random coefficients and a minimal parametrization of the compact 
Stiefel manifold related to these random coefficients. Using this parametrization of the 
general class, a complete identification procedure has been proposed. Such a statistical 
inverse problem requires to solve the stochastic boundary value problem in high stochas- 
tic dimension with efficient and accurate algorithms. New results of the mathematical 
and numerical analyzes of the parameterized stochastic elliptic boundary value problem 
have been presented. The numerical solution provides an explicit approximation of the 
application that maps the parameterized general class of random fields to the correspond- 
ing set of random solutions. Since the proposed general class of random fields possibly 
contain random fields which are not uniformly bounded, a particular mathematical anal- 
ysis has been developed and dedicated approximation methods have been introduced. 
In order to obtain an efficient algorithm for constructing the approximation of this very 
high-dimensional map, we have introduced complexity reduction methods using low-rank 
approximation methods that exploit the tensor structure of the solution which results 
from the parametrization of the general class of random fields. 
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